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ABSTRACT 

We present the analysis of two Suzaku observations of GX 301-2 at two orbital phases after the periastron 
passage. Variations in the column density of the line-of-sight absorber are observed, consistent with accretion 
from a clumpy wind. In addition to a CRSF, multiple fluorescence emission lines were detected in both ob- 
servations. The variations in the pulse profiles and the CRSF throughout the pulse phase have a signature of a 
magnetic dipole field. Using a simple dipole model we calculated the expected magnetic field values for differ- 
ent pulse phases and were able to extract a set of geometrical angles, loosely constraining the dipole geometry 
in the neutron star. From the variation of the CRSF width and energy, we found a geometrical solution for the 
dipole, making the inclination consistent with previously published values. 

Subject headings: X-rays: stars — X-rays: binaries — stars: pulsars: individual (GX 301-2) — stars: magnetic 
fields 


1. INTRODUCTION 

The High Mass X-ray Binary (HMXB) system GX 301-2 
was discovered in 1969 April during a balloon experiment 
(Lewin et al. 1971; McCIintock, Ricker & Lewin 1971). The 
system consists of an accreting neutron star (NS) fed by the 
surrounding stellar wind of the B type emission line compan- 
ion Wray 977 (Jones, Chetin & Liller 1974). A recent lu- 
minosity estimate derived from atmospheric models puts its 
distance at ~ 3 kpc (Kaper, van der Meer & Najarro 2006), 
the value utilized in this paper. The orbital period was es- 
tablished to be ~ 41 days (White, Mason & Sanford 1978) 
using Ariel 5 observations and was refined with the Burst 
And Transient Source Experiment (BATSE) to ~ 41.5 days 
with an eccentricity of ~ 0.46 (Koh et al. 1997). Doroshenko 
et al. (2010a) discussed a possible orbital evolution and de- 
termined an orbital period of 41 .482 ± 0.001 d, assuming no 
change in orbital period. Kaper, van der Meer & Najarro 
(2006) determined that the mass of the companion was in the 
range 39 Mq < M < 53 Mg and the radius of Wray 977 was 
R* ~62R q , obtained by fitting atmosphere models. 

The X-ray flux is highly variable throughout an individ- 
ual binary orbit but follows a distinct pattern when averaged 
over multiple orbits (see Figure 1). Shortly before the peri- 
astron passage, the X-ray luminosity increases drastically in 
the energy band above ~ 5 keV, as seen in Rossi X-ray Tim- 
ing Explorer (RXTE)/ All Sky Monitor (ASM) data (Leahy 
2002). The NS passes closest to the companion at a distance 
of ~ 0. 1R» (Pravdo et al. 1995). Shortly after the periastron 
passage, <p 0 * ~ 0.2, the X-ray luminosity dips for a short pe- 
riod of time. Leahy (2002) demonstrated that neither a simple 
spherical wind model nor a circumstellar disk model around 
Wray 977 are sufficient to describe the observed variations in 
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FIG. 1.— RXTE/ASM 1 .5- 1 2 keV light curve folded with the orbital pe- 
riod of 41 .48 days. The two lines indicate the times of the orbital phase when 
the Suzaku observations were performed. The profile shows the strong flux 
increase during the pre-periastrion flare at orbital phase ~ 0.9. For clarity, 
two binary orbits are shown. 

the folded RXTEIASM data. An additional stream compo- 
nent is able to account for the sudden increase in X-ray lumi- 
nosity, as the NS passes trough the stream shortly before peri- 
astron and accretes more material. This model also explains a 
slightly higher X-ray luminosity around o 0r b ~ 0.5, when the 
NS passes through the accretion stream a second time. 

Pulsations with a period of ~ 700 s were discovered in the 
Ariel-5 observations (White et al. 1976), making GX 301-2 
one of the slowest known pulsars. The pulse period has varied 
drastically throughout the last ~ 20 years (Pravdo & Ghosh 
2001 ; Evangelista et al. 2010). Prior to 1984, the pulse period 
stayed relatively constant at 695 s-700 s and then spun up be- 
tween 1 985 and 1990 to ~ 675 s. From 1 993 until the begin- 
ning of 2008, the change in the spin reversed again, showing a 
decline. Fermi/Galactic Burst Monitor (GBM) data 6 have re- 
vealed that GX 301-2 experienced another spin reversal and 

6 http://www.batse.msfc.nasa.gov/gbm/science/ 




2 


Suchy et al . 


briefly spun-up with the pulse period decreasing from ~ 687 s 
to ~ 681 s between May 2008 and October 2010. Since Oc- 
tober 2010, the pulse period has shown only small variations 
around ~ 681 s. 

Most recently, Gogtij, Kreykenbohm & Belloni (201 1) dis- 
covered a peculiar 1 ks dip in the luminosity of GX 301-2, 
where the pulsations disappeared for one spin cycle during 
the dip. Several such dips have been previously observed in 
Vela X-l (Kreykenbohm et al. 1999, 2008), where it is as- 
sumed that the accretion on the NS was interrupted for a short 
period of time. 

The pulse phase average spectrum of GX301-2 is de- 
scribed using a power law with a high energy cutoff. The 
continuum does not show a strong variation in the intrinsic 
parameters (r,£ cut , and £f 0 id) throughout the orbit (Mukher- 
jee & Paul 2004), as seen in two data sets from RXTE, taken 
in 1996 and 2000, sampling most phases of the binary or- 
bit. One of the major characteristics of the X-ray spectrum 
of GX 301-2 is the high and strongly variable column density 
of its line-of-sight absorber throughout the orbit, indicative 
of a clumpy stellar wind (Ah = 10 22 - 10 24 cm 2 ). In addi- 
tion to the high column density, a very bright Fe Ka emission 
line can be observed. This line has shown a strong corre- 
lation with the observed luminosity, indicating that the line 
is produced by local clumpy matter surrounding the neutron 
star (Mukherjee & Paul 2004). Kreykenbohm et al. (2004) 
used the RXTE data set from 2000 to perform phase resolved 
spectroscopy and showed that an absorbed and partially cov- 
ered pulsar continuum (power law with Fermi-Dirac cutoff) 
as well as a reflected and absorbed pulsar continuum were 
consistent with the data. 

A cyclotron resonance scattering feature (CRSF) at ~ 
35keV was first discovered with Ginga (Mihara 1995). Or- 
landini et al . (2000) found sy stematic deviations from a power 
law continuum at ~ 20 and ~ 40 keV in BeppoSAX , where 
the former could not be confirmed as a CRSF due to the prox- 
imity of the continuum cutoff. Kreykenbohm et al. (2004) ex- 
cluded the existence of a CRSF at ~ 20 keV and showed that 
the CRSF centroid energy varies between 30-38 keV over the 
pulse rotation of the NS. Furthermore, they showed that the 
CRSF centroid energy and width are correlated. 

We report on two observations of GX 301-2 performed 
with the Suzaku satellite in mid 2008 and early 2009. This 
paper is structured as follows: Section 2 discusses the obser- 
vations and data reduction. Section 3 shows the phase aver- 
aged results. Section 4 discusses the pulse profiles and phase 
resolved spectra. Sections 5 and 6 discuss the results and con- 
clusions, respectively. 

2. OBSERVATION AND DATA REDUCTION 

Suzaku observed GX 301-2 on 2008 August 25 with an ex- 
posure time of ~ 10 ks (ObsID 403044010; hereafter Obs. 1 ). 
The observation was cut short by a set of target of opportunity 
observations and was continued on 2009 January 5, acquir- 
ing an additional ~ 60 ks exposure time (ObsID 403044020; 
Obs. 2). Both main instruments, the X-ray Imaging Spectrom- 
eter (XIS; Mitsuda et al. 2007) and the Hard X-ray Detector 
(HXD; Takahashi et al. 2007) were used in these observations. 
The two observations correspond to orbital phases of 0.19 
(Obs. 1) and 0.38 (Obs. 2), where Obs. 1 falls into the lowest 
flux part of the binary orbit (see Figure 1). The 2- lOkeV 
absorbed flux was 1 .6 x 1 CT 10 erg cnrT 2 s _1 for Obs. 1 and 
~ 8 x 1 0~ 10 erg cm " 2 s _l for Obs. 2 (see Table 1 for details). 
Both observations were performed using the HXD nominal 


pointing to enhance the sensitivity of the HXD detectors. 

The XIS detectors consist of two front illuminated (FI) 
CCDs XIS 0 and 3, and one back illuminated (BI) CCD, XIS 
1 . All three instruments are sensitive between ~ 0.5 keV and 
~ 1 1 keV, although the two FI cameras have a higher sensi- 
tivity above ~ 2keV, and the BI chip is more sensitive below 
~ 2keV. To minimize possible pile-up, the XIS instruments 
were operated with the 1/4 window option with a readout time 
of 2 s. Data were taken in both 3x3 and 5x5 editing modes, 
which were extracted individually with the Suzaku FTOOLS 
version 16 as part of HEASOFT 6.9. The unfiltered XIS data 
were reprocessed with the most recent calibration files avail- 
able and then screened with the standard selection criteria as 
described by the Suzaku ABC guide 7 . The response matrices 
(RMFs) and effective areas (ARFs) were weighted according 
to the exposure times of the different editing modes. The XIS 
data were then grouped with the number of channels per en- 
ergy bin corresponding to the half width half maximum of the 
spectral resolution, i.e. grouped by 8, 12, 14, 16, 18, 20, and 
22 channels starting at 0.5, 1 , 2, 3, 4, 5, 6, and 7 keV, respec- 
tively (M. A. Nowak, 20 10, private communication). The XIS 
spectral data were used in the energy range of 2- lOkeV for 
reasons described below. 

The HXD consists of two non-imaging instruments: the 
PIN silicon diodes (sensitive between 12-60keV) and the 
GSO/BGO phoswich counters (sensitive above ~ 30 keV). To 
determine the PIN background, the Suzaku HXD team pro- 
vides the tuned PIN Non X-ray Background (NXB) for each 
individual observation. In addition, the cosmic X-ray Back- 
ground (CXB) was simulated following the example in the 
ABC guide and both backgrounds were added together. The 
PIN data were grouped by a factor of 5 below and by a fac- 
tor of 10 above 50keV for Obs. 1. In Obs. 2 the channels 
were grouped by 3 throughout the whole energy range. GSO 
data were extracted and binned following the Suzaku ABC 
guide. The PIN data energy range of 15-60keV and the GSO 
data energy range of 50-90 keV were used for the phase aver- 
aged analysis. Due to the short exposure times, no GSO data 
were extracted in Obs. 1 and in the phase resolved analysis of 
Obs. 2. 


2.1. XIS Responses 

GX 301-2 is a source with very large and highly variable 
photoelectric absorption in the line of sight as well as a very 
strong Fe Ka emission line. For spectra] modeling, the cal- 
culated responses of the XIS detectors showed a ‘leakage’ 
emission level (Matsumoto et al. 2006), which stemmed from 
the instrument characteristics. That is, a small fraction of an 
event charge cloud is registered at energies below the peak en- 
ergy, creating a ’low energy tail’ in the spectrum. This ‘tail’ 
is about three orders of magnitude weaker than the intensity 
of the peak count rate of the measured photon. The response 
to the Fe line is therefore described as a Gaussian plus a con- 
stant at lower energies, approximately three orders of magni- 
tude below the peak value and extending to below 1 keV. Due 
to the strong Fe Ka line and the very strong absorption in 
this source, the constant level response of the Fe line extends 
above the measured flux below 2 keV. This results in signifi- 
cant residuals below ~ 2 keV, which seem to be stronger for 
the BI spectrum then for that of the FI CCD. Figure 2 shows 
the phase averaged spectrum from the second observation of 

’ http:/7heasarc .gsfc .nasa.gov/docs/suzaku/anaIysis/abc/ 
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FIG. 2.— (a) Best fit models of the Suzaku GX301-2 second observation 
for the 3 individual XIS instruments, and the two HXD instruments. At lower 
energies, one clearly sees the constant level in the modeled flux discussed in 
the text, (b) shows that the residuals below ~ 2keV are very pronounced. 

GX 301-2 and the best fit models using the HEASOFT gen- 
erated XIS responses as detailed in section 3. 

Together with the Suzaku Guest Observer Facility (GOF), 
which provided for us an experimental response parameter file 
where the 'low energy tail’ is ignored, we studied this behav- 
ior in more detail. This new response for all three instruments 
was applied to the individual XIS spectra. Due to the fact 
that the new response completely excluded the missing "tail’, 
we now observed an excess of the data with respect to the 
model below ~ 2keV. The introduction of additional power 
law and/or black body components could improve the fits at 
lower energies. Due to the fact that the true level of the tail is 
not known, these additional components cannot be interpreted 
physically, and we decided to use the original response matrix 
and confine the selected XIS energy range to above 2 keV. 

3. PHASE AVERAGED SPECTRUM 
3.1. Spectral modeling 

Phase averaged broad band spectra in the 2-60 keV (Obs. 1) 
and 2-90 keV (Obs. 2) energy ranges were obtained. From the 
technical description of the XIS instruments 8 , it is known that 
small discrepancies, e.g., in fitted power law slope, have been 
observed between FI and BI XIS instruments. A difference 
of ~ 0.05 in the power law index has been observed in cali- 
bration data and was also observed here. When modeling the 
FI and BI instruments with a common power law index, the 
residuals of the BI instrument show a systematic deviation. 
However, only the FI or BI instruments could be modeled to- 
gether with the HXD instruments. Due to this fact and the 
higher sensitivity of the FI XIS instruments above 2keV, we 
concentrated our discussion on the results obtained with the 
two FI instruments. Each data set, FI and BI, was modeled 
individually with the HXD data to compare the differences in 
the continuum. We found that when using the FI and BI in- 
struments individually with the HXD instruments, in both ob- 
servations the best fit spectral values are consistent with each 
other within error bars, indicating that the usage of the HXD 
reduced the observed discrepancies in the power law parame- 

8 http://heasarc.nasa.gov/docs/suzaku/prop_tools/suzaku_td/ 


ters. 

The continuum model for pulsars can so far only be mod- 
eled with empirical models, consisting of a powerlaw with 
a cutoff at higher energies, which is typical for this kind of 
source (Cobum 2001). Three empirical models are widely 
used: the simple high energy cutoff (highecut) and the 
Fermi-Dirac cutoff (fdcut, Tanaka 1986) both are used in 
combination with a simple powerlaw component. In addition, 
the negative-positive exponential powerlaw model NPEX (Mi- 
hara 1995) is a slightly more complicated model including the 
power law component. For the data analyzed in this work, the 
best results have been obtained with the fdcut model: 

E~ r 

/corn— A pl — — 7= 7 , (1) 

exp(£-E cut /£f 0 id)+l 

where Apt is the normalization at 1 keV, £ cut is the cutoff en- 
ergy, and £foid is the folding energy of the Fermi-Dirac cut- 
off. For the Obs. 1 and Obs. 2 FI data sets the best fit x 2 ’s 
were 260 and 542 with 245 and 267 degrees of freedom (dof), 
respectively (see Figure 3 and 4). Replacing the smoother 
Fermi-Dirac cutoff with the highecut model in Obs. 2,re- 
sulted in slightly worse residuals (A^ 2 « 1 5) compared to the 
best fit results of the FD-cutoff model (Figure 3 d) and an 
emission line-like residual at ~ 35 keV, which is most likely 
due to the sharp break of the power law at the cutoff energy 
in this model. In comparison, the NPEX model shows even 
worse residuals (Ax 2 ~ 200) and only results in reasonable 
fits when the exponential curvature at higher energies is inde- 
pendent in the partially and fully covered component. 

The best fit cutoff energy is very close to that of the ob- 
served CRSF feature (see below) and these two parameters 
are rather strongly correlated (see Figure 5). To avoid a de- 
generacy of the CRSF values due to a changing cutoff energy, 
E C ut was frozen in all the fits to the best fit value from the 
Obs. 2 FI spectrum (29.2 keV). 

Previous observations of GX 301-2 indicated the existence 
of clumps in the stellar wind (Kreykenbohm et al. 2004; 
Mukherjee & Paul 2004), which were modeled using partial 
covering absorption in addition to fully covered photoelectric 
absorption of the smooth stellar wind. In the present analysis, 
the low energy portion of the spectrum also required a partial 
covering component (Figure 3b and 4b), which was modeled 
using the TBnew model (Wilms et al., 2011, in prep.) 9 , an 
updated version of the existing TBabs model (Wilms, Allen 
& McCray 2000). In addition, a non-relativistic, optically- 
thin Compton scattering component cabs was included, as 
is necessary for column densities ;Vh > 5 x 10 22 cm' 2 , where 
the plasma becomes Compton thick and part of the emission 
is scattered out of the line of sight. The Ah values from 
TBabs and from the cabs component were set equal and 
treated as one model component for the smooth stellar wind 
(Ah.i = TBabs l*cabsl) and a second component for the 
clumpy partial coverer (A H ,2 = TBabs2*cabs2). The Ah.i 
and Ah ,2 column densities were left independent of each other. 
For spectral fitting the abundances of Wilms, Allen & McCray 
(2000) and the cross-sections (Verner et al. 1996) were used 
in all data sets. 

Both observations show residuals to the continuum in the 
35-40keV energy range, which we interpret as the pre- 
viously observed CRSF (Mihara 1995; Kreykenbohm et al. 
2004). Modeling these residuals with an absorption line with 
a Gaussian optical depth (gabs) improved the residuals sig- 

9 http://pulsar.stemwarte.uni-erlangen.de/wilms/research/tbabs/ 
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Fl G . 3 . — Data and best fit model for Obs. 1 . a) shows the data and the 
best fit model. The inset shows the residual for the different emission lines 
observed in XIS data, where the Fe Ka line is already included in the model, 
b) shows the best fit residuals without partial covering and the CRSF. c) has 
only the CRSF added, and d) has both componenst included. 

nificantly in both observations (see Figures 3 and 4). The best 
fit values of £crsf = 42. 5 keV and 35. 2^; 9 keV for Obs. 1 
and Obs. 2, respectively, are consistent within 90% confi- 
dence intervals. The widths of the CRSFs, ocrsf = 9.5*|;f keV 
(Obs. 1) and 7-8Zq A 9 keV (Obs. 2), are also consistent within 
errors. An increase of ctcrsf with higher centroid energy, 
as previously observed in RXTE data (Kreykenbohm et al. 
2004), could not be detected, although the best fit values hint 
at such a behavior. Table I shows the best fit values for the 
continuum with a gabs CRSF line where the given errors are 
90% confidence values. To calculate the significance of the 
CRSF in the first and weaker observation, the null hypothe- 
sis approach was applied, where 10,000 spectra were created 
with Monte Carlo simulations using the best fit parameters 
without the CRSF. Gaussian uncertainties were used for the 
individual model parameter. Each model was fitted with and 
without the CRSF line to compare how much an inclusion of 
the line actually improves the fit. In ~ 99.6% of all fits, no 
CRSF feature were observed with a larger optical depth than 
the observed lower limits of r = 10.3, making the existence 
of the feature significant with a ~ 3. Including the CRSF in 
the real data, the best fit improved by a A\ 2 = 1 00. A similar 
improvement of \ 2 was only observed in ~ 0.1 % (a ~ 3.3) of 
the simulated spectra, concluding that the observed line in the 
first observation is indeed real. For the second observation, 
the line was even more pronounced, where the exclusion of 
the CRSF increased the x 2 by ~ 200. 

The CRSF can alternatively be described with the Loren- 


FlG. 4.— Same as Figure 3 for Obs. 2. In this case b) shows the residuals 
without partial covering and the CRSF. c) has only the CRSF added and d) 
has both components included. 
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Fig. 5.— Contour plots of £crsf vs Cut ■ The different contours indicate 
the 1,2, and 3 sigma contours for the FI spectrum of Obs. 2. The X indicates 
the best fit values for both parameters: Ecu = 29.2 keV and fcRSF = 35.2 keV. 

zian shaped eye labs XSPEC model (Mihara et al. 1990). 
Cyclabs is described by the centroid energy £crsf, the 
width er CRSF , and the resonance depth tcrsf. similar to the 
gabs parameters. The best fit continuum parameters are 
consistent with the best fit values determined with the gabs 
component. The cutoff energy with a best fit value of 
£■0111 = 31 .2^'9keV is consistent with the value obtained with 
the gabs model. The observed centroid energies of £crsf = 
35.0^2jkeV and 31.6^;°keV are of the order of 10-20% 
lower than the energies obtained with the gabs model. This 
discrepancy stems from a different calculation of the line cen- 
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TABLE 1 

PHASE AVERAGED SPECTRAL PARAMETERS with the gabs model component for best fit spectrum 


Parameter 

403044010 (lOks) 

403044020 (60ks) 


FI 

BI 

FI 

BI 

Wh.iIIO 22 cm” 2 ] 

16.62f? 

16.6 (frozen) 

20.9+ 3 f 

20.9 (frozen) 

Abund Ca 

1 .55 (frozen) 

1 .55 (frozen) 

1 cc+O.39 
1 - > ~’-0. 18 

1 9g+0.41 
1 - yo -0.49 

Abund Fe 

1.17 (frozen) 

1.17 (frozen) 

1 t 7 +o.04 
1,1 7 -0.03 

1 04+0.03 
1 ,U ^-O.04 

yVjf- 110 22 cm” 2 ] 

76.9 1 *” 

76.9 (frozen) 

28.4!j-“ 

28.4 (frozen) 

r 

n q^+o.08 

U.o3_o g£ 


0.965;“; 

0 9C+0- 03 
U - VO -0.03 

Apm [ i 0“ 3 j 1 fully covered 

1 68 40 - 67 
1 • oo -n.50 

1 -30+0.20 

1 ■'*9-0.20 

1 1 +u 
1 * 1 -0.5 

4 7 o+°- 9 i 
/u -0.79 

Apl ,2 [Hr 3 ] 1 part, covered 

92 S + 16 - 1 

90.02| 6 , 2 

215 . 55 - 9 1 

266.42!” 

Ecut [keV] 

29.2 (frozen) 

29 2 (frozen) 

29.2 + I 3 ~ > 

29.2 (frozen) 

Efold |keV| 

10-3!^ 

8.3f 2 “ 

5.7®? 

5 . 6 !?;} 

IScrsf [keV] 

42.5«J 

39.7 +11 

35-2 

34.921? 

ccrsf [keV] 

o 5 + 1.6 
- 2.8 

8.65-0 

1 r 41 - 1 

7 4 + 0-5 
' *^-0.5 

tcrsf 

29 3 +17 - 7 

18.9f‘ 4 / 

12 . 0! 22 , 7 

10.62}-* 

£skc [keV] 

9 oo+O.Ol 

7 71+0.03 
1 -0.03 

7 - 37 +0.01 
+'^_ 0.01 

7 34 +O.bl 
^^- 0.01 

IskoJIO ” 5 ] 2 

9 -14+0.61 

? 44+0.96 
^•^-0.89 

5 36+0.50 

3 " 3O -0.50 

4 94 +O .68 
^•^-0.67 

EQWj Ka [eV] 

860^ 24 

1 109 +1313 
1 1U *-947 

671^1 

965^ 2 

Ea rKa [keV] 

9 qa;+ 0-03 
^• yt> -0.03 

3 02 4<) ()6 

‘ > - u +'_0.05 

3 OO 4 ^ 1 
• 3 - UU - 0.02 

3.°4*j;<p 

4 53+0.94 

/ArKcJlO - 5 ] 2 

i 79 +O .62 
1 ' '^-0.69 

1 OCT 0 - 94 
1 * uu -0.48 

2.74!«i 

EQW 3 rKct [«V] 

179!^ 2 

1 3 0 +3()8 
1 1U -1 10 

38tf® 

642 3 * 

£ca Ka [keV ] 

-3 70+0.05 

J * /J -0.04 

3-« 

3 70 + ^ > 02 
-QQ 1 

3 71+0.03 

/caKallO ” 5 ] 2 


2-45!};i 

4-72t^ 

4 95 +I .68 
-1.71 

EQW| aKQ leV] 


l0 °i!m 

14^8° 

!4!!i 

EpcKa [keV] 

6 399 + 0.007 
°" )W 4(X)6 

6 410 +()-()l0 

6 409 40 003 

6+24 

O-FcKa [eV] 

13! 22 

<31 

<5 

<7 

/ftkallO ” 4 ] 2 

7.16!®“ 


‘9.20!S;i 

99 9+0.61 
^ ,v -0.58 

EQWpcKa [eV] 

241 2 g 

230 4 ^ 

13323 7 

1422}? 

EfeK0 [keVJ 

7 06+ 003 
' • UO -0.03 

7 Od 4 ^- 94 
' - W -0.04 

7.0920;?,! 

7 09+0- 0 - 
' ,uy - 0.02 

/FeKtfPO - 4 ] 2 

1 19 + 0.29 
1 • 1 *-035 

1 54 + 0.99 
1 •~* + -0.60 

2.65^;g 

4 5Q+0.53 
^•-^-0.56 
32+9 
^ -10 

EQWjUa [eV] 

432“ 

56+® 

212 2 

£ni Ka [keV] 

- 

- 

7 46 4 ^ 02 
'•^-o.oz 

' ' J 1 -0.06 

hi Ka 1 10 ” 4 ] 2 

- 

- 

1 57 + 0 . 2-1 
1-3/ -0.25 

°-672?:l 

EQW 3 iKQ [eV] 

- 

- 

152? 

5.72.i 3 7 4 

Rux f lOteV absorb ' 

I-6H8S 

1 63 +0 - (B 
1 - DO -0.07 

7.98®-°? 

9 O4+0.04 
y,U ^-0.06 

Flux^_ 10keV unabsorb. 

16.3! 2 ;? 

16 O 408 

1 O - U -0.7 

31.82“? 

OC Q+l .() 

’ 53 - y -0.4 

Cxis 3 /C P ,N/Cg S0 

0.93^/1.18!«;«/- 

-/ 1 12 40 - 07 / - 
' 1 ' 1 A-0.06 ' 

0 95 l' (4)003 /i 31 + 0 . 02/1 39 + 0.12 
* _o .003 / 101 -0.02/ 1 -' ? ^-0.12 

_ / 1 1 7 + 0.02 / 1 9 7+0- 1 1 
‘ ia / -0.02/ 1 * Z/ -0.11 

X 2 /dof 

260 / 245 

145/ 130 

542 / 267 

328/ ISO 


(l)Units are ph keV” 1 cm 2 s (2) Units are ph cm 2 s 1 , (3) Values of EQWs are determined relative to the abs. continuum, 


(4) Absorbed and unabsorbed flux units are [10 10 erg sec 1 cm 2 ], (5) Values of C are with respect to XIS 0 for the FI fits and XIS 1 for the BI fits. 


troid energy and is described in detail in Nakajima, Mihara 
& Makishima (2010). The width ctcrsf = 1 1.22§QkeV and 
12.31} 3 keV forObs. 1 and Obs. 2, respectively, is bigger than 
with the gabs model. The x 2 I dof values of 260 / 245 and 
548 / 267 for the two observations using the cyclabs model 
and thus were not a better fit when compared to the gabs 
model. For the final discussion we use the values determined 
by the gabs model. 

Several emission features were observed in the residuals of 
the fits in both observations and were subsequentially mod- 
eled with Gaussian emission lines (see Figures 3 and 4, inlay). 
The width of each line was set equal to that of the Fe Ka emis- 
sion line, while the intensities were left to vary independently. 
Energies were loosely constrained around the expected values 
of neutral material to avoid runaway of the line energies. 

A constant (const) was applied, taking small instrumental 
differences in the overall flux normalization into account. The 
constant was fixed at 1 for XIS 0 and was left free for the other 
instruments. C\is 3 , Grin and Cgso in Table 1 are the cross 
calibration constants with respect to XIS 0 for FI fits and with 
respect to XIS 1 for the BI fits. 

The final model had the form: const*fVH,i (PL1+ 


N Hi 2 *PL 2 )* f dcut*gabs+ 6 xGaussians, where the 
power law indices of PL1 and PL2 are set to be equal to each 
other. The individual power law normalizations were inde- 
pendent and were used to calculate the fraction of the partial 
covering. 

Table 1 summarizes the best fit values for both observations 
and for the individual FI / BI data sets. For the interpretation 
we will concentrate on the FI data for both observations, for 
reasons mentioned above. 

3.2. Spectral results 

The best fit spectral parameters were generally consistent 
with previous observations, with the exception of the cutoff 
energy. With a value of 29.221 3 ; 5 keV, E ctA was significantly 
higher than the ~ 20keV value measured in previous RXTE 
(Mukherjee & Paul 2004) and BeppoSAX (La Barbera et al. 
2005) observations, although these observations used slightly 
different spectral models for the cutoff energy. Kreykenbohm 
et al. (2004) used also the here applied Fermi-Dirac cutoff, 
resulting in cutoff energies of 10 - 15 keV. The folding energy 
for Obs. 1. £f 0 id = 10.3ij QkeV, is consistent with values ob- 
tained for a similar orbital phase with BeppoSAX (La Barbera 
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et al. 2005). In Obs. 2, the value of Zifoid = 5.7 ^ 5 keV is con- 
sistent with the RXTE data for the pre-periastron flare and the 
periastron passage (Kreykenbohm et al. 2004). 

The observed column densities for the absorption due to 
the smooth stellar wind (Mhj) were consistent between both 
observations. A larger difference could be observed in the 
column density associated with the absorption due to the 
clumped wind (Afu), where the best fit value of Obs. 1 is 
almost a factor three larger than that of Obs. 2. In addition 
to the H column density in the TBnew model, the relative 
Ca and Fe abundances were also left independent in Obs. 2. 
The best fit values are measured from the Ca and Fe K edges 
at 4.1 keV and 7.1 keV, respectively, and were slightly higher 
than solar abundances: 1 .55^| for Ca and 1 -17^j ^ for Fe. 
For the first observation, the Ca and Fe abundances could not 
be well constrained and were fixed to 1 .55 for Ca and 1 . 17 for 
Fe. 

The partial covering fractions can be calculated from the 
measured normalization values of the two power law compo- 
nents, Apu andApL 2 : 

Cvr. Frac. = — — . (2) 

ApLl +ApL 2 

The factor is almost unity in both observations due to the dom- 
inance of the Apl 2 value. Similar large covering fractions have 
been observed in XMM-Newton data (Fiirst et al. 201 lb) and 
in the RXTE data (Mukherjee & Paul 2004), indicating that 
the covering fraction does not change significantly for differ- 
ent parts of the orbit. The spectrum softened slightly from 
T = 0.83^ for Obs. 1 to 0.96!°$ for Obs. 2. 

The Fe Ka line was detected at 6.409^^ keV (Obs. 2), 
which is, given the instrumental gain systematics, consistent 
with neutral Fe. In addition, the following lines have also 
been observed in both observations: S Ka: line, Ar Ka: line, 
Ca Ka line, Fe K j3 line, and the Ni Ka line, all with energies 
consistent with neutral material (e.g. Kaastra & Mewe 1993). 
Note that the Ni Ka line was not detected in the fainter first 
observation. The observed line intensities are summarized in 
table 1 with 90% confidence errors. The values are consistent 
between the FI and BI instruments. The widths of the lines 
were set to be equal to the Fe Ka width within each observa- 
tion, which was found to be 13l| 2 eV for Obs. 1 and had an 
upper limit of 5 eV in Obs. 2 (FI values). The Compton shoul- 
der to the Fe Ka line, as observed by the Chandra observatory 
(Watanabe et al. 2003), was not significantly detected in either 
observation. The observation of Watanabe et al. (2003) was 
performed in the pre-periastron phase, where the luminosity 
was higher than in the two Suzaku observations. A Compton 
shoulder has also been observed with XMM-Newton in data 
taken during another pre-periastron flare (Fiirst et al. 201 lb). 

4. PHASE RESOLVED ANALYSIS 

Barycentric and binary corrected light curves for differ- 
ent energy bands were extracted for both observations us- 
ing the orbital parameters from Koh et al. (1997) with an 
updated orbital period and periastron time t 0 determined by 
Doroshenko et al. (2010a). Due to the strong variability 
of the pulse period on short time scales, individually deter- 
mined pulse periods were used for each observation. An ac- 
curate pulse period could be established for Obs. 2, but not 
for Obs. 1 due to its short duration of lOks. For the sec- 
ond observation the calculated pulse period had a value of 
P = 685 .4 ± 0.9 s. The Fermi/GBM instrument measures the 
pulse period of GX 301-2 on a regular basis, resulting in val- 
ues of P = 687.28 s for Obs. 1 and P = 685.75 s for Obs. 2 



Pulse Phase 


Fig. 6.— Pulse profiles for different energy bands for Obs. 2. The two 
peaks are indicated in the top panel . The bottom pannel shows the selection 
of phase bins which were used for phase resolved spectroscopy, 

(Finger. 201 1 , priv. comm). To create the pulse profiles, the 
Fermi IG B M-prov ided pulse periods and epoch times for each 
observation were used. 

We studied the energy dependent pulse profiles for Obs. 1 
and Obs. 2. Both pulse profiles are very similar, showing a 
double peak shape where the main peak (PI) is broader below 
lOkeV. The second peak (1*2) stays rather constant in width, 
but increases its relative intensity toward higher energies. As 
an example, the pulse profile for Obs. 2 in different energy 
bands is shown in Figure 6. Comparing these pulse profiles 
with previous RXTE and BeppoSAX data showed that the 
general shape is consistent through all parts of the orbit. 

The data of Ihe longer Obs. 2 were divided into 10 equally 
spaced phase bins and individual spectra were extracted for 
the XIS and PIN instruments. Figure 6 shows the individual 
phase bins in the lowest panel. This division resulted in an 
exposure time in each phase bin of ~ 5 - 6 ks for each individ- 
ual instrument. No GSO data were used in this analysis due 
to the reduced exposure time per phase bin. For spectral anal- 
ysis, the phase averaged model was applied to all phase bins, 
resulting in best fit values summarized in Table 2, as well as 
Figure 7. Again the cutoff energy was frozen to the phase 
averaged value of 29.2 keV to avoid the previously discussed 
degeneracy with £crsf- 

The two absorbing components, Af H .i from the smooth stel- 
lar wind and ;V H ,2 from the partial coverer vary between 20 
and 40 x 10 22 cm -2 , but show an anti-correlated trend through- 
out the first peak. Although the error bars are rather big, one 
can see an indication that A/h,i follows the flux in PI , whereas 
the (Vh ,2 value dips at the same time. 

In contrast to Kreykenbohm et al. (2004), the power law 
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TABLE 2 

Phase resolved spectral parameter fromObs. 2 for FI data. Same model as phase averaged data. 


Parameter 

PBl 

PB2 

PB3 

PB4 

PBS 

PB 6 

PB7 

PBS 

PB9 

PB10 

AklllO 22 /™! 2 ] 

19.2«°7 2 

20 . 0 !™ 

36.12}® 

70 ? +23 - 6 

32-45f 3 

22.8 2™ 

25.32*4 

12 . 223 - 

71 4+ 7 -9 
•‘+- 7.3 

24.62™ 

W n , 2 it 0 22 /cm 2 I 

32.62™ 

30.92+] 

71 3+13-8 

7 5 n+ 21 - C * 
- 20.6 

18-75i 

26.72™ 

32.7+1] 

35.2+20 

29 5 +3 - 9 

-3.9 

31.6211 

r 

1 - 22 3 * 

1 30+°^ -fr* 

1 OU -0.04 

0 93+0-04 
u * VJ> -0.04 

0 91+0-03 
U,yi -0.04 

1 7 Q+ 0 - 0 4 

1 *“ U -0.04 

1.07!«-‘0 

0 74 + 0.05 
U - 7 M).05 

0 56 +<) - 04 
u - 3 °-o m 

o-675:g 

1 1 c+0.06 
1 • 1 ' J -0.06 

ApL, 1 [ 1 0 “ 2 ] 1 ful . cov . 

°-^:™ 

1 35 + 3 - 13 

10. 125 2 } 

0-295f 

8.69*fg 7 

1.922};5 

2 03 +1 53 
- 1.00 

0 . 162 O-® 

o-®5:£ 

2.362}'®} 

ApL^l lO ” 2 ] 1 par. cov. 

24.42}] 

44.8^1 

1 7 4 + 6.6 

27.4!t 9 

41. D™ 

27.751 

10.82™ 

10 - 621:11 

1 1 .3+14 

19.325'} 

Gold fkeVJ 

7.4!‘V 

6 13+ 0 - 38 
0 - 1 -5-0.35 

5-6^ 

5 S * 0 5 

D ' D -0 A 

59 + 0. 7 

3 -^- 0.6 

4 1+5.0 


5 0 + 0.5 

d - u -oa 

4 4 + 0.3 
M).3 

5 0 + 0.4 
3 - U -0.3 

Ecrsf [keV] 

38.7*-” 

29.8}jJ 

32.lt™ 

34.654 

38.CF|] 

38.3 5™ 

35-4 2 > : > 

35.7+jO 

31 1+ 44 

97 5 +I.I 

<rcRSF [keV] 

10.35} 

5 r 1 - 7 

6-4 5;? 

7.354 

6.3 !|; 2 

15 . 15 ™ 

8.42 jO 

8-i5:? 

6-75:! 

4 9 + 1 . 1 
^•‘"- 0.8 

TCRSF 

33 4+75.1 
-16.7 

3 5 + 1.5 

6.71 2 j 

8.45] 

8 . 35-4 

30 . 52^:8 

17.0 + ‘ o A 

21-855 

n.42 2 J 

0 q+ 1.2 
J * ~-0.9 

feKaflO - 3 ] 2 

i Q7+0.09 

1,y/ -0.09 

1 87+0.10 
L * Z -0.10 

1 -^-o.os 

1 89^- 13 
1 ’ oy -0.07 

1 79 + 0.09 
1 * /Z -().08 

i-762S:1t 

i-87t!f 

7 16 40 ' 11 
^• lo - 0 . 1 l 

2 ll * 0 - 51 
z * 11 -o.u 

7 07+0- 10 
‘~ u/ - 0.10 

RuX tl0keV 

A Q3+0.03 
" K7 ~ , -0.08 

7 77 +O.O 8 
'-0.16 

1 0 03+ 4 - 42 
1U - UJ> -0.61 

1 1 97+0.39 
,Zr -3.37 

1 0 91+ 4 37 
iU -^°-G.90 

7.762 

6-565:1 

§ 47 + 0.17 
o>f/ - 0.20 

7 A 7 -K). 10 
' • Uj “- 0 . 2'2 

4 89 f0X,:> 

CxlS 3 /Cp Is . 

0.94/1.15 

0.94/1.15 

0.95/1.19 

0 . 94 / 1.21 

0 . 94 / 1.21 

0.94/ 1.19 

0.95/1.16 

0.95/ 1.14 

0 . 94 / 1.10 

0 . 94 / 1.11 

X 2 /dof 

233 / 243 

288 / 243 

353 / 243 

312/243 

253 /243 

303 / 243 

283 / 243 

292 / 243 

264 / 243 

247 / 243 

X 2 /dof no CRSF 

241 / 246 

307 / 246 

391 / 246 

341 / 246 

273 / 246 

305 / 246 

291/246 

344 / 246 

297 / 246 

270 / 246 


(1) Units are ph keV 1 cm 2 s 1 , (2) Units are ph cm 2 s ',(3) Units are 10 10 erg cm 2 sec 1 , (4) Values of C are with respect to XIS 0 for the FI fits. 


index T is varying strongly throughout the pulse. The val- 
ues are 0.9-1 .3 throughout the first peak and drop suddenly to 
0.6— 0.8 for the second peak (Figure 7). Measuring the second 
peak to be significantly harder than the first peak supports the 
behavior seen in the pulse profile, where the intensity of P2 in- 
creases at higher energies, becoming similar to the intensity in 
PI . Power law normalizations are relatively small and badly 
constrained for the first power law component. The partially 
covered power law normalization seems to follow the flux, 
showing a higher value throughout the first peak. The folding 
energy does not vary significantly throughout the orbit and 
shows values between 5 and 6 keV (not shown in Figure 7) 

The CRSF energy varies between 30-+10 keV, similar to 
the values observed in the RXTE data (Kreykenbohm et al. 
2004). Table 2 shows that the addition of the CRSF in the in- 
dividual spectra improved the x 2 for most of the phase bins, 
except in phase bin 6, which falls in the gap between both PI 
and P2. Phase bins 1 and 7, also minima in the pulse profile, 
only show a small improvement in the best fit when the CRSF 
was included. The CRSF energy is not very well constrained 
in phase bins I and 6 due to the lack of sufficient statistics. 
For the geometrical discussion below, these two phase bins 
are ignored. In all other phase bins the addition of the CRSF 
improves the fit. The energy changes very smoothly through- 
out the pulse and does not follow explicitly the observed flux 
in each phase bin. The CRSF energy increases during the 
main peak and reaches the maximum at the falling flank of 
PI . During P2 the CRSF energy decreases until it reaches a 
minimum at the dip between P2 and PI . 

To estimate the significance of the CRSF in the phase re- 
solved analysis, again the null hypothesis method was applied 
in two out of the 8 remaining phase bins, which are good rep- 
resentatives of all data points. Phase bin 2 is a good example 
for a shallow CRSF, whereas phase bin 5 is an example for a 
phase bin with a higher flux. Similarly as with the phase aver- 
aged data, 10,000 spectra were simulated and modeled with- 
out and with the CRSF. The interpretation of the observed 
feature as being due to stochastic fluctuations with a depth 
similar to the observed values could be ruled out with a prob- 
ability of over 99% (a ~ 3) in both phase bins. This leads to 
the conclusion that the CRSF can be observed with sufficient 
significance to allow the interpretation below. 

5. DISCUSSION 


5.1. Phase averaged continuum 

In both observations very strong absorption can be ob- 
served. A?h,i , which is the absorption due to the smooth stellar 
wind, is consistent within both observations. Ah. 2 , the absorp- 
tion from the partial coverer, is significantly stronger in the 
first observation, establishing the existence of an additional 
component in the line of sight. Leahy & Kostka (2008) used 
RXTE / ASM and PCA data to show that a stellar wind and a 
stream model component can describe the observed count rate 
throughout the orbit and that the increase of the column den- 
sity at orbital phase ~ 0.2 is mainly due to the stream compo- 
nent. The findings in both Suzaku observations are supportive 
of this theoretical picture. 

GX 301-2 shows very strong luminosity variability on very 
short time scales (Kreykenbohm et al. 2004; Ftirst et al. 
2011b), where the column density increases to up to 50% in 
periods of lower activity. Both RXTE and XMM-Newton 
data showed a much higher column density in the pre- 
periastron flare than in the two data points observed with 
Suzaku. One possible explanation is that the wind is in- 
deed so variable and clumpy that these variations are just not 
properly predictable, especially in times of high activity, such 
as the pre-periastron flare. Leahy & Kostka (2008) investi- 
gated RXTE / PCA Ah values of the absorbed component from 
archival observations for different parts of the orbit and found 
no increased column density in the pre-periastron flare. Ob- 
servations throughout one full binary orbit could help to un- 
derstand the variations in Ah on time scales of days. 

The observed luminosity dependence of the power law in- 
dex F and the folding energy is similar to other sources, such 
as V0332+53 (Mowlavi et al. 2006) and 4U01 15+63 (Tsy- 
gankov et al. 2007). The power law index T shows a hard- 
ening for the first observation, whereas the folding energy 
£foi d is slightly higher when compared to Obs. 2. Soong et al. 
(1990) observed the £f 0 ] d variation in Her X-l phase resolved 
HEAO-1 data and concluded that the parameter is dependent 
on the viewing angle of the accretion column and can be used 
to describe the plasma temperature of the system. In the case 
of GX 301-2 the smaller £f 0 id value in Obs. 2 could indicate 
a lower plasma temperature which can be interpreted that the 
X-ray emission region is further up in the accretion column 
(Basko & Sunyaev 1976) (see also Section 5.2). 

The softening of the power law index T with increased lu- 
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Pulse Phase 

Fig. 7.— Phase resolved spectral parameters for Obs. 2. Top panel shows 
the 20 -40 keV pulse profile, where the individual phase bins are numbered. 

Units of power law normalization are ph keV 1 cm" 2 s" 1 . The intensity of 

the Fe line is units of ph cm" 2 s" 1 . For clarity, two pulses are shown. 


minosity is also in agreement with the basic model of the ac- 
cretion column that the plasma temperature decreases with in- 
creased height. With increased luminosity, the rate of accreted 
material would increase and the amount of soft photons cre- 
ated by the lateral walls of the relatively taller column would 
increase, leading to a softer spectrum, as is observed. 

For the second observation the abundances for Fe and Ca 
were left independent and a slight overabundance was ob- 
served. Taking into consideration that the abundances used 
in Wilms, Allen & McCray (2000) are derived from the inter- 
stellar medium (ISM) a small overabundance from an evolved 
star may be expected . 

5.2. Variations in the CRSF parameters 

In both observations the CRSF was clearly observed and 
the gabs absorption component improved the overall fit sig- 
nificantly. Strong magnetic fields (B ~ 10 12 G) exist close to 
the NS magnetic poles, where photons at energies close to the 
Landau levels are resonantly scattered from the line of sight 
and result in an absorption line-like feature in the spectrum. 
This feature provides a direct method to measure the magnetic 
field strength close to the NS surface, where the fundamental 
centroid energy can be described as 

Bcrsf = 1 1 .6 keV x x (3) 

where z is the gravitational redshift at the scattering site, with 
z « 0.3 for typical NS. Using the Ecrsf values obtained in the 
phase averaged spectra, the calculated magnetic field has a 
strength of 4.76^ x 10 12 G(Obs. 1) and 4.10!^ x 10 12 G 
(Obs. 2) which is consistent within errors for both observa- 
tions. 

The CRSF parameters are consistent with previous results 
using RXTE data and are lower than the observed BeppoSAX 
(La Barbera et al. 2005) and INTEGRAL (Doroshenko et al. 
2010b) values of 45-50keV. These observations were ob- 
tained during pre-periastron outburst, where the luminosity of 
the source is much higher than in the two Suzaku observa- 
tions presented in this paper. A luminosity dependence of the 
phase-averaged CRSF centroid energy is observed in multi- 
ple other sources, where an anti-correlation between CRSF 
and luminosity was observed in V 0332+53 and 4U01 15+63 
(Tsygankov et al. 2006; Mihara et al. 2007), and a posi- 
tive correlation was observed in Her X-l (Staubert et al. 
2007). Klochkov et al. (2011) have confirmed, from pulse 
to pulse variability studies, such correlations for V 0332+53, 
4U01 15+63 and Her X-l, and also found an anti-correlation 
for A 0535+26. Although the CRSF values arc consistent 
with a single value in both Suzaku observations, a hint of 
an anti-correlation can be seen in the data. Compared with 
the BeppoSAX and INTEGRAL data, the CRSF centroid en- 
ergy might have a positive correlation. Parallel to the CRSF- 
luminosity correlation, Klochkov et al. (201 1) also observed 
in pulse to pulse data of four sources that the power law index 
F shows an opposite luminosity dependence than the Bcrsf- 
In the Suzaku observations discussed here, F is observed to 
soften with increased luminosity, another indication that a 
possible negative CRSF-luminosity correlation exists. 

The two opposite scenarios seem to depend on the luminos- 
ity, i.e., if the observed luminosity is below or above a critical 
luminosity (CL) which can be derived from the Eddington lu- 
minosity of a system (Becker & Wolff 2007 ). The CL is of 
order ~ 1 0 j7 erg s" 1 , but depends also on the accretion geom- 
etry (spherical or disc) as well as the height and diameter of 
the accretion column. 
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Above the CL, the infalling proton density becomes so large 
that the protons begin to interact and decelerate, creating a ra- 
diation pressure dominated shock region above the magnetic 
pole. This region of increased density is most likely the region 
where the CRSF is created (Basko & Sunyaev 1 973). With in- 
creasing luminosity, the shock region moves higher up in the 
accretion column, where a smaller local magnetic field value 
results in a lower observed CRSF centroid energy, as is ob- 
served in V 0332+53 or 4U01 15+63. The observed CRSF- 
luminosity dependence as well as the T-luminosity depen- 
dence in GX 301-2 is very similar to the scenario described 
here. 

Below the CL, the accreting matter slows down via hy- 
drodynamical shock, 'Coulomb friction’ or nuclear collision 
(Basko & Sunyaev 1973; Braun & Yahel 1984). With in- 
creasing accretion rate, the luminosity increases and the de- 
celeration region is pushed closer to the NS surface, where 
the magnetic fields are higher. This would result in a posi- 
tive CRSF-luminosity correlation, as observed in Her X-l and 
A 0535+26. 

With a distance of 3 kpc , the intrinsic unabsorbed 2-10 keV 
luminosity of ~ 2 x 10 33 ergs _1 is significantly below the typ- 
ical CL of ~ 10 37 ergs _l . Although it has been observed that 
sources with similar luminosities can show the opposite cor- 
relation (Klochkov et al. 201 1) due to individual critical lumi- 
nosity values, the large luminosity difference compared to CL 
puts this source in the same regime as Her X- 1 and A 0535+26 
and at odds with the notion of a well-defined CL. The ob- 
served variation in T and possibly in £crsf would then stem 
from another mechanism. 

5.3. Emission lines 

The existence of multiple fluorescence emission lines, es- 
pecially at lower energies, where the absorption is very domi- 
nant, indicates that the source of the emission originates from 
a region where the column density is not very large, i.e. the 
outer layers of the stellar wind in our line of sight (Fiirst et al. 
201 lb). If the emission lines would be embedded deeper in 
the stellar wind, the lines, especially at lower energies would 
have to be significantly stronger, to be detected at all. On the 
other hand, if the line emitting region is in the outer layers 
of the wind, the incident soft X-ray flux would be drastically 
reduced, making the equivalent width very large, as observed 
in the low energy emission lines (Table 1). A possible expla- 
nation would be that the emission region is very large for the 
lines, and maybe spread over the entire surface of the stellar 
wind. 

The most dominant line emission stems from the Fe Ka 
transition at 6.4 keV. The observed energies in both observa- 
tions are consistent with the emission in neutral material. The 
ratio of the intensities of Fe K 3/Fe Ka for both observations, 
0.1638* for Obs. 1 and 0. 143;®j for Obs. 2, is consistent with 
the fluorescing material being neutral or only slightly ionized 
(Kaastra & Mewe 1993). The equivalent widths (EQW) of the 
two Fe Ka lines show that the intensity relative to the contin- 
uum is reduced by a factor of ~ 2 in the second observation. 
Based on emission line widths, Endo et al. (2002) estimated 
in ASCA data that the Fe emission originates within 10 lo cm 
of the continuum emission source. Fiirst et al. (201 lb) use 
XMM-Newton pulse by pulse observations to observe a cor- 
relation with increased flux of the source, showing that the 
Fe emission region is not far from the X-ray source. In the 
Suzaku observation, however, the Fe Ka line flux does not 
directly follow the observed luminosity. The increase in ab- 


sorbed 2-10 keV luminosity is a factor of ~ 5 , whereas the 
increase in the fluorescence line intensities is only 2-3. This 
difference in intensity change indicates that the distance be- 
tween continuum and fluorescence emission lines is large, es- 
pecially when compared to the proposed ~ 0.3 lt-s from Endo 
et al. (2002). In the 6 months’ time between the two obser- 
vations, however, the overall emission geometry could have 
changed, accounting for the difference in the fluorescence in- 
tensities. If the emission originates from a confined region 
in the stellar wind, observations taken during two different 
orbital phases would each yield a different distance to this re- 
gion. 

In addition to the strong Fe Ka line, many other emission 
lines are observed in the XIS spectra in both observations. 
The observed intensities of these lines show the same behav- 
ior, where the flux for the second observation is ~ 50% or 
less compared to the first observation, although the observed 
uncertainties for the different lines, especially for S Ka and 
Ar Ka, are very large. Furthermore, the S Ka- line falls into an 
energy regime where calibration uncertainties have been pre- 
viously observed (e.g., Suchy et al. 201 1), so that a detailed 
analysis is not possible at this time. 

5.4. Pulse profiles 

Historically, GX 301-2 shows very strong pulse to pulse 
variations and an intensive flaring behavior on very short time 
scales (Kreykenbohm et al. 2004; Fiirst et al. 2011b), e.g., 
throughout the pre-periastron flare. The observations obtained 
here do not show such flaring behavior, and only show the 
regular pulsation throughout the observation. A pulse period 
could only be established for Obs. 2 of the Suzaku data, which 
was consistent with Fermi/GBM data. With the GBM pulse 
period for the first observation, pulse profiles for both obser- 
vations could be produced for different energy bands. The 
two peaked pulse profiles show a similar behavior for both ob- 
servations, where the second peak is getting strong at higher 
energies. This behavior is in contrast to many similar sources, 
such as 4U 0115+63 (Tsygankov et al. 2007), V 0332+53 
(T sy gankov et al . 2006) , 4U 1 909+07 (Fiirst et al . 20 1 1 a) , and 
1 Al 1 18-61 (Suchy et al. 201 1), where the pulse profile turns 
into a single peak profile at higher energies. A reduction in 
the intensity of the second pulse has also been observed in 
BeppoSAX data for different parts of the orbit (La Barbera 
et al. 2005). 

5.5. Geometrical constraints using a simple dipole model 

The phase resolved analysis of the second observation 
showed strong variations in several spectral parameters 
throughout the pulse profile. As previously observed by 
Kreykenbohm et al. (2004), the CRSF centroid energy varies 
by ~ 30% where the highest energy is detected at the peak 
and the falling flank of PI . Such behavior has been observed 
in multiple other sources, e.g., Vela X-l (La Barbera et al. 
2003), 4U 0352+309 (Cobum 2001), and Her X-l (Klochkov 
et al. 2008). A model for the CRSF variations is based upon 
the change in the viewing angle throughout the pulse and thus 
different heights of the accretion column are probed, yielding 
a different local magnetic field observed for each phase bin. 

A simple approach to derive a possible geometry for the 
neutron star and the magnetic field uses the variation in the 
observed magnetic field throughout the pulse phase. In this 
case the very smooth and sinusoidal variation (Figure 8) can 
be modeled as a simple dipole where the total magnetic mo- 
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FIG. 8.— Calculated phase resolved magnetic fields from the CRSF energy 
and the best fit for a simple magnetic dipole model. 

ment jj is calculated from the phase averaged CRSF energy 
value. 

The variation of the CRSF energy is then fitted by chang- 
ing the geometrical angles of the system until a best fit con- 
verges. Appendix A discusses the model in detail , introducing 
the three free parameters: 0, the viewing angle between line 
of sight and the NS spin axis, the inclination angle a of the 
magnetic moment with respect to the spin axis, and the an- 
gle 3, indicating the ’lag’ of the B-field plane with respect 
to the ephemeris, i.e., the observed shift in pulse phase (see 
also Figure 9) The best fit values for ,3 were either -4.6° or 
175.4°, depending on the starting values of the fit (see Ap- 
pendix A), where the latter value is equivalent to the former 
when the dipole polarity is flipped, a and 0 angles showed a 
very strong interdependence, and the best fit 0/a pairs were: 
-14.6° /67°, 14.6°/113°, 67°/- 14.6°, and 113°/14.6°. All 
best fit values have a x 2 ~ 4 with 7 dof. These pairs do show a 
degeneracy of 0± 14.6° and 90±23°, which can be explained 
with a geometrical symmetry, when calculating the magnetic 
field for each individual phase bin. The first and last two pairs 
can be each treated as the same geometry, rotated by 180 de- 
grees. Figure 8 shows the calculated magnetic field values 
from the CRSF centroid energies and the theoretical values 
for B from the best fit angles for a simple dipole model with 
0 = 67°, a =-14.6° and ,3 = -4.6°. 

The phase resolved energy and width of the CRSF show a 
strong correlation, where the width varies in phase with the 
CRSF energy (Figure. 7). Kreykenbohm et al. (2004) ob- 
served a similar correlation, where the magnetic pole was ob- 
served under different viewing angles <I> , where <t> is the angle 
between the line of sight and the magnetic axis. Meszaros 
& Nagel (1985) showed that the anisotropic velocity field of 
the electrons in the accretion column leads to a fractional line 
width of: 

fe . 1 *-? . 1 . ~ \l 81n2— — r jcos^| =£jcos$|. (4) 

£crsf V m e c 2 

The <tcrsf/£crsf ra hos vary in the 0.16 - 0.24 range, cor- 
responding to variations of ±20% throughout the pulse phase. 
The two outliers with large errors in phase bins 1 and 6 were 
ignored. The variation of 4> can be estimated as 0 -q and 
0 + a throughout one pulse phase, where the magnetic pole 
rotates around the NS spin axis, tilted by 0 with respect to 
the line of sight. 

Using the average ratio of 0.2, we calculated the values of k 
for the two geometries with 0 = -14.6° (k = 0.2 1) and 0 = 67° 
(k = 0.51). For each geometry, we used the variation of the 


ocrsf/£crsf ratio (±0.04) to calculate the variation in 4>, re- 
sulting in values of $ = 15° ± 25° and <f> = 67° ± 5°. These 
variations are smaller then the values expected from the the 
geometrical discussion, but show a similar behavior, where 
the variation of the angle is larger for the smaller value of <J>. 
From Equation 4, we were also able to estimate a possible 
plasma temperature for the two geometries from the calcu- 
lated values of k. The plasma temperature of kT ~ 4keV for 
0 = 15° is very similar to the observed folding energy fifoid, 
which is an indication of the plasma temperature (Burderi 
et al. 2000, and references therein). The value of kT ~ 24keV 
for the 0 = 67° geometry is much larger and is not consistent 
with the observed £f 0 id- Assuming that the NS spin axis is 
aligned with the inclination of the binary system, a value of 
0 = 15° would put the inclination at i ~ 75°, a value which is 
only marginally above the upper limit of 72°, as determined 
by Kaper, van der Meer & Najarro (2006). 

The strong Fe line was detected in all 10 phase bins. The 
measured line flux did not change significantly throughout the 
pulse phase, while the 2- lOkeV flux varied by more than a 
factor of 2. From this we conclude that the distance to the Fe 
fluorescence region is greater than ~ 700 lt-s (~ 2 x 1 0 13 cm), 
from the NS. This is in conflict with the conclusion of Endo 
et al. (2002) that the emission region must be closer than 
10 lo cm from the NS, based upon the width of the emis- 
sion lines measured by ASCA. Their assumption is that the 
emission region is close to the NS and that the material free 
falls onto the NS surface, encountering much faster veloci- 
ties, which caused the observed broadening of the lines. The 
smaller measured Fe Ka line width in the phase averaged 
spectra is more consistent with the assumption that the line 
broadening is due to the terminal velocity of 300 -400 km of 
the line driven wind (Parkes et al. 1980). 

6. SUMMARY 

We presented results of two Suzaku observations of 
GX301-2 taken shortly after the periastron passage. The 
spectra were modeled with a partially and fully covered power 
law with a Fermi-Dirac cutoff and a CRSF at ~ 35keV. 
The column density of the smooth stellar component only 
marginally changed and the clumpy wind component signifi- 
cantly changed between both observations. Flux dependen- 
cies of the CRSF and T were not significant, but do hint 
at similar correlations such as those observed in V 0332+53 
and 4U01 15+63, although the observed luminosity is signifi- 
cantly below the calculated critical luminosity of this system. 
They are in the realm where one might expect the energy of 
the CRSF to decrease with decreasing flux. The variations in 
the pulse profiles and the CRSF throughout the pulse phase 
have a signature of a dipole magnetic field. Using a simple 
dipole model we calculated the expected magnetic field val- 
ues for different pulse phases and were able to extract a set of 
geometrical angles, loosely constraining the dipole geometry 
in the NS. Model constraints derived using the observed ratio 
of ^ crsf to £crsf, together with calculated plasma tempera- 
tures, favor a solution wherein the spin axis is tilted ~ 15° to 
our line of sight. 
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FIG . A 1 . — Schematic view of the dipole model, indicating the different angles used to calculate the CRSF energy variation .0 is the angle between the projection 
of the magnetic moment in the x-y-plane and the ephemeris (x-axis) .The eye indicates the line of sight towards us. 

Suzaku satellite, a collaborative mission between the space agencies of Japan (JAXA) and the USA (NASA). 

APPENDIX 

A. GEOMETRICAL DIPOLE MODEL 

The variability of the CRSF energy throughout the pulse phase can be used to understand the geometry of the neutron star. We 
assume a simple dipole magnetic field around a spherical neutron star, however we want to emphasize that these assumptions are 
very rudimentary, and effects such as the displacement of the dipole field from the neutron star center are not taken into account 
here. 

The total magnetic moment /r of the neutron star can be calculated with: 

(n) = B ■ Rj5, s ( vT+3cos^¥) = 6.776 x 10 3O G cm 3 (Al) 

using the average magnetic field from the phase averaged CRSF centroid energy (B = 4. 1 x 10 12 G) and a typical neutron star 
radius of R^s = 10 km. The value of the last part, including the azimuthal angle <f> of the magnetic pole was estimated to the 
average value of 1 .54. 

The total average magnetic moment is divided into three components based on the cartesian coordinate system in Figure Al: 

( sin(a)cos(,3)\ 

sin(a)sin(/3) 1 (A2) 

sin(;3) J 
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where a is the angle between the dipole field and the spin axis and / 3 is the rotation angle from the ephemeris, i.e. the x-axis. 
Next, the vector n was calculated for each individual phase bin (!) used in the phase resolved analysis. 


( sin(0)cos(d>)\ 
sin(0)sin(<p) | 
sin(< 2 >) ) 


(A3) 


The individual components of n are the projections of each phase bin on the three axes, as seen from the line of sight of the 
observer. The angle 0 corresponds to the angle between the line of sight and the spin axis. For each of the individual phase bins, 
the magnetic field was calculated using the position vector n indicating the position of the magnetic pole and the vector m, which 
indicates the magnetic moment: 


3-n|m-n| m 


(A4) 


The calculated magnetic field values were then fitted to the measured magnetic field, varying the three angles 0, a, and (3 until 
the best fit converged. To avoid wrong solutions due to local minima, the starting parameter of all three angles were varied in 
steps of 10° before the fit was started. The best fit solutions resulted in the values that are discussed in detail in section 5.5. 
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